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ABSTRACT 


A new capability has been added to the general purpose finite element 
program NASTRAN Level 17.7 to conduct forced vibration analysis of tuned cyclic 
structures rotating about their axis of synrretry. The effects of Coriolis and 
centripetal accelerations together with tliose due to linear acceleration of the 
axis of rotation have been included. 

This report presents the theoretical development of this new capability. 
The work was conducted under Contract NAS3-22533 from NASA Lewis Research 
Center, Cleveland, Ohio, with Mr. Richard E. Morris as the Technical Monitor. 
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SUMMARY 

The objective of the work described herein, was the development, docu- 
mentation, demonstration and delivery of a computer program for the forced 
vibration analysis of rotating cyclic structures. Tuned bladed discs are an 
example of specific interest. 

The scope V'^as by definition to address: 

- direct periodic loads moving with the rotating structure , 
specified in the frequency or time domain; 

- translational acceleration of the rotating axis. 

The capability is operational in the NASTRAM general purpose program at 
Level 17.7. 

flASTRAN documentation is provided and example analysis results have been 
obtained. 

Relationships to previous work are described and further developments 
are reconmended. 
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1 . INTRODUCTION 


In Reference 1 (NASA CRT 59728), a general approach to conducting dynamics 
analysis of bladcd discs was discussed and a logical sequence of problems defined. 
Figure 1 indicates the problem spectrum end Figure 2 shows an overall program 
structure of modular nature, based on additions and modifications to the NASTRAN 
general purpose structural analysis program. 

The general set of problems can be represented by formal equations, various 
terms of which are included or excluded depending on the problem being studied. 

Thus the equations 

[MJ {ii} + [[B] + 2n[B^H{u) - iuZ (u) + [[Kj - [H,]] {ii} 

= (P) - [M2n’Ro} > (1) 

with some boundary conditions, may be taken to represent the general forced 
vibration dynamic problem of a tuned bladed disc. In the previous and current 
work, the advantages of cyclic symmetry are incorporated in this formulation. 

The work of Reference 1 handled, within state-of-the-art technology, analyses 
of aeroelastic, modal and flutter problems of tuned cyclic structures (Figure 2). 

It was documented and operated in NASTRAN at Levels 16 (Refs. 2 - 5), and 17.7 (Ref. 6). 

The equations 


[MJ {ijl + [BJ {u} - [QJ (u) + [K] (u) = {0} (2) 

were treated in the context of modal and flutter problems of a tuned system, 
where 

H represented the inertia matrix; 

B represented the damping matrix; 

Q represented the induced aerodynamic matrix (complex); 
and K represented the elastic cum geometric stiffness matrix. 

In the current work, the equations 


[HJ {u} + [[ 6 ] + 2'a [B^J]{u} + [[K] - [H^J] {u} = {P) - {Rq} 

(3) 


are treated in the context of forced vibration, (Figure 2) > 
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Figure 2; Overall Program Structure & Status 
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where [Ml, [B], [Kl are as before, 

[B-|] represents Coriolis forces matrix 

■[Pj represents appl ied surface load vector 

[M2l represents base forcing mass matrix 
and {Rq} represents base acceleration vector. 

The specific form of the forced vibration equations including the types of 
forcing functions to be incorporated were selected by NASA. These are specifically 

- directly applied loads moving with tne rotating structure; 

- inertial loads due to translational accelerations of the ax’s of 
rotation ("base acceleration"). 

The leads may be periodic and specified in the frequency or time domains. 
Solution procedures follow generally those of the cyclic symmetry formulation 
of the NASTRAN Theoretical Manual (Ref. 7). The capability has been develo ad 
on the IBM 370 system at Bell Aerospace Textron and documented, delivered and 
demonstrated on the UNIVAC 1108 at NASA Lewis in NASTRAN Level 17.7. 

Five demonstration examples are presented, one being a complete structure 
example to shew compatibility with the cyclic structure formulation. A simple 
twelve bladed disc is modelled and forced at conditions related to its natural 
frequencies. The response examples include; 

- physical component forcing (frequency domain), 

- harmonic component forcing (frequency domain), 

- harmonic component forcing of the rotational axis, 

- physical periodic forcing (time aomain), 

- harmonic component periodic forcing (time domain) 

It would be logical to extend the current v.ork to include the generation of 
applied and induced oscillatory aerodynamic loads sc hat the forced vibrations 
of subcritical engine stages can be addressed directly. Rotational base accel- 
erations could also be of practical interest in investigating the gyroscopic 
effects on rotating machinery. 
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2. EOli'ATICHS OF MOTION 


The equations of motion of a tuned cyclic structure rotating about its axis 
of syirsTietry and subjected to steady sinusoidal and general periodic excitation 
are derived using Lagrange's formulation. 

Figure 3 illustrates the problem by considering a 12-bladed disc as an 
example. The bladed disc consists of twelve 30° segments — identical 

in their geometric, material and constraint properties. The :.isc rotates about 
its axis of symmetry at a constant angular velocity. The axis of rotation itself 
is peniiitted to oscillate linearly in any given inertial reference. In addition, 
the bladed disc is allowed to be loaded with steady sinusoidal or general periodic 
loads moving with the structure. Under these conditions, it is desired to deter- 
mine the dynamic response (displacement, acceleration, stress, etc.) of the 
bladed disc. 

The cyclic symmetry feature of the rotating structure is utilized in deriving 
and solving the equations of forced motion. Consequently, only one of the cyclic 
sectors is modelled and analyzed using finite elements, yielding substantial 
savings in the analysis cost. Results, however, are obtained for the entire 
structure. The Coriolis and centripetal acceleration terms have been included. 

For clarity of derivation, the equations of motion are first derived for an arbi- 
trary grid point of the cyclic sector finite element model, and then extended 
for the complete model. 

COORDINATE SYSTEMS 

These are shown in Figure 3. 0-XYZ is an inertial coordinate syste i. 
0-XgYgZg is a body-fixed coordinate system such that OXg coincides with the axis 
of rotation of the structure and is always parallel to OX. For a NASTRAN finite 
element model of the bladed d’sc, 0-XgYgZg also represents the ^asic coordinate 
system. A-xyz is a body-fixed global coordinate system in which the displace- 
ments of any grid point P are desired. The unit vectors associated with these 
coordinate systems are also shown in Figure 3 . 

DE GREES OF FREED OM 

The rotating structu>^e is permitted four rigid body motions including three 
translations (along OX, 6Y and OZ) and one rotation at a constant angular velocity 
n about its axis of rotation OXg. 



y 



I, 3 , K Unit vectors along Inertial XYZ axes 

^'E vectors along Basic XgYgZg axes 

i, j» k Unit vectors along Global xyz axes 


Figure 3: Coordinate Systems 
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OF POOR QUALITY 


All grid points of the structure are each permitted six degrees of freedom. 
The displacement at any grid point in any sector can be expressed in any body- 
fixed coordinate system as a combination of: 

1) the steady state displacement due to the steady rotation of, and the 
steady state loads applied to, the structure, and 

2) the vibratory displacement (superposed on the steady displacement) due 
to the vibratory excitation provided by the directly applied loads and the 
inertial loads due to the acceleration of tiie axis of rotation ('base' 
acceleration). 

The purpose of the present development is to determine the vibratory 
response. 

LAGRANGE FORMULATION 

Referring to Figure 3, the complete tuned structure consists of N identi- 
cal cyclic sectors. If u represents all the vibratory degrees of freedom of the 
complete structure, the equations of motion can be derived via the Lagrange 
formulation. 


dt ' aa ' ■ 3u 3u 38 ■ 3u 


where T and U represent the kinetic and strain energies, respectively, of the 
complete structure; D is the Rayleigh's dissipation function representing the 
energy lest in the system due to resisting forces proportional to velocities 3 
(e.g. viscous damping forces); and 6W represents the virtual v;ork done on the 
structure by the external forces through virtual displacements 6u. 


The complete set of degrees of freedom u can be subdivided into N subsets, 
each containing u'^ degrees of freedom for each of the N cyclic sectors. Since 
any given cyclic sector is 'connected' to adjacent cyclic sectors only on its 


two sides, u" satisfies the intersector boundary compatibility condition 


u 


n 

side 2 


u 


,n+l 
side 1 


n = 1, 2. 


N. 


{?.) 


Equations (1), therefore, can be written as N sets of equations coupled 
only as given by equations (2): 


^ , ;)T" . 3T" , 3U" . ^ 

dt ' 3un ^ ^ aflfi " ’ 


n = 1, 2, 


(3) 
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For clarity of presentatiorj, without loss of generality, equations (3) are 
first applied to obtain the equations of motion of an arbitrary grid point in any 
cyclic sector by considering its three translational degrees of freedom. Inclu- 
sion of the three rotational degrees of freedom at the arbitrary grid point, and 
extension to include the remaining grid points in the cyclic sector are con- 
sidered subsequently. 

KINETIC ENERGY 

4* h 

With reference to Figure 3, point P is an arbitrary grid point of the n*^" 
cyclic sector with a mass of 'm' units lumped from the adjacent finite elements. 


The kinetic energy of the mass at P can be written as: 


where 


and 


with 


T = 2 m Pp • rp 


^^P ^OP ^ ^ ^OP^’ 

o O /s o ,s o 

^0 ■ * V * 


^ O ^ O O 

PQP "" "OP^B ^OP'^B * ^0P*^B ’ 


= ni , 


(>] 

3d > 


B 

A\ 


1 0 0 

0 c s 

0 -s c 


c H cos nt and s e sin Jit 


Substitution of equations (5) through (9) in equation (C) results in 
1 


0 0 0 


T = i LXq Yq ZqJ 


'm 

/ o \ 

( ^o) 


\ 0 / 

m 

Yo> 


/ 0 \ 

m 

V^oj 

— — 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


(9) 


(10) 



' m 


o \ 

*0P 

O O 0 

L^OP ''^OP ^opj 

m 


> 


m 


V^op j 
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In order to introduce the global coordinates of point P, consider now the position 
vector to P written as 


0? = OA + A? , (12) 
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Therefore equation (13) yields 






OF POOR QUALITY 


noting that (16) 

The global position vector to the point P can be further thought of as consisting 
of three. component vectors, i.e.. 


+ \ u.. > + ( u 


AP 

initial 


Initial pos ition a Pis 
bisplacements due 
to n and steady 
loads 


Displacements cue to vibratory forces 


Steady part 


Vibratory part 


Substitution of equations (17) in equations (15), and differentiation once with 
respect to t yields 


[tGB^ 


1 

u_ 


Therefore, with the help of equations (15), (17) and (18), equation (11) expresses 
the kinetic energy of the mass at P in terms of the displacements and velocities 
of P expressed in the global coordinate system A-xyz. 

STRAIN ENERGY 

The strain energy due to the displacements at point P expressed in the basic 
coordinate system can be written as 

U = ^ [_(Xqp - Xgp^ , (Yqp - Yqp^ .^_) , (7qp - Zqp^ * 


*^XX "^XY "^XZ 


'YX YY 


L ^ZX ^ZY •''ZZ 


^YZ ) ' 
K- 7-7 Z 


^OP ■ ^OP.in. -j 

Y - Y ( 

^OP ’op, in. ^ 

V . z J 

*-0? '^OP.in. 
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By using equations (15) and (17), this becomes 


op FCOfJ Q'.JAUTf 


Global 


l!., + U 


1 I — — — uluyai i _ 




U H- U J 
Z Z 


where 


Global -n T PD 

[K J = . 


( 20 ) 


( 21 ) 


DISSIPATED ENERGY 


The energy dissipated by the damping forces acting against the motion of the 

point P is expressed by the Rayleigh's dissipation function, 

Bisic 




Ynn Znn_l 


2 L.'‘op 'OP '-0P 


Bu„ B 


Use of equation (18) transforms this to 


1 I o o o I Global 

^ = “zJ ^ 


u 

<" U ") 

I;: 

'• II 


B. 


XX “XY “XZ 

Byx Byy By2 

®ZX , ®ZY ^ZZ 


OP 

OP 


‘^OP 


(22) 


(23) 


where 


Global PD T Basic po 
[B J = irh [BJ irh 


(24) 


VIRTUAL WORK 

The virtual work done by an oscillatory applied load 
^ = ^x’ V ’’z^ 


(25) 


moving with the point of appi ication, through individually achieved virtual 

displacements 6u , 6u and 6u can be stated as 
X y 2 


6W 


5 LiJx % Pv 


(26) 
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EQUATIONS OF FORCEO i TjTIOH 

Substitution of the expressions for T (equation 11), U (equation 20), D 
Equation 23) and 6W (equation 26) in the Lagrange equations (3) results in the 
following equations of forced motion of point P expressed in the displacement 
(global) coordinate system A-xyz: 


[M] {ii} + 


[B] -f- 2fi[B^J 


(u) + 


The terms appearing in equations (27) are: 


[K] - (u) = (P) - [M^J {Rjj} 


(27) 


{u} = { u 


'-1 


'■ U. 


(28) 


r P, 

(pi'l p. 


L 


(29) 


{Rq} = Yq 


(30) 


( r 
P I 


global -n T 
[HJ = [M J = [T^^J 


global 
[BJ = [B 2 


0 

in 

0 


0 

0 

m 


[tGBj 


(31) 


(24) 


global 

[B^J = [B^ ] 


global 
[KJ = [k J 


[t“j 


0 0 
0 -m 
m 0 




(32) 

( 21 ) 


global 


[H^J = [M^ J = [T^^J 


0 

m 

0 


0 

0 

m 


^tGBj 


(33) 
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end 


[M2l 


y ni 0 0 

[T^^] 0 me ms 

0 -ms me 


> 


where c and s are given by equation (10). 


(34) 


Equations(27) describe the translatory motion of an arbitrary point P in 
an arbitrary sector n of the rotating cyclic structure subjected to a directly 
applied vibratory load (P) and base acceleration (Rq). 

These equations can be extended to include the three rotational degrees of 
freedom at point P by noting that: 

1) in a lumped mass model, only the translational degrees of freedom at 
any grid point contribute to the kinetic energy of the structure, and 

2) the coupling between various degrees of freedom may exist only via 
the stiffness matrix. (Instances where the damping matrix is defined 
proportional to the stiffness matrix also may result in coupled equations 
of motion. ) 


viz. 

by 


Accordingly, the matrices derived from kinetic energy considerations, 
[Mj, [B-j], [M^J and [M 2 ] of equation (27) can be expanded as typified 


[M] 

6x6 


I _ 0 

0 '1 0 


(35) 


where [M^j_] is the 3x3 (^translational ) mass matrix of equation (27). With sub- 
scripts t and r representing the Jtranslational and Rotational degrees of freedom 
at point P, the stiffness and damping matrices may be expanded as 


(36) 


and 


[K] = 

1 

't-I 

1 

1 


.f^tJ i 

CK,rJj 


■[Btt] ; 

CBtrl ■ 

II 

1 — I 

CO 


* — ' ~ 



[E,rJ 


(37) 


By similar reasoning, the equations of forced vibratory motion of all the 
cyclic sectors of the total structure can be written as 
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The intersegment boundary compatibility is specified by equation (2). 
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3. SOLUTION OF EQUATIONS OF MOTION 


The method of solution of the equations of forced motion (equations 38 and 2 , 
Section 2 ) is based upon the form in v/hich the excitation of the rotating struc- 
ture is specified. As noted earlier, the present development considers excitation 
prescribed as: 

1) directly applied loads moving with the sti'ucture and 

2) inertial loads due to the translational acceleration of the axis of 
rotation ('base' acceleration). 

These steady-state sinusoidal or general periodic loads are specified to 
represent: 

1) the physical loads on various segments of the complete structure, cr 

2) the circumferential harmonic components of the loads in (1). 

The sinusoidal loads are specified as functions of frequency and the general 
periodic loads are specified as functions of time. 

The translational acceleration of the axis of rotation is specified as a 
function of frequency in an inertial coordinate system. 

Because of its eventual implementation in the HASTRAN general purpose finite 
element structural analysis program, the following solution procedure generally 
follows the theoretical presentation of cyclic symmetry given in the NASTRAN 
Theoretical Manual (Ref. 7). 

METHOD OF SOLUTION 

The method of solution of the equations of motion consists of four principal 
steps: 

1) Transformation of applied loads to frequency-dependent circumferential 
harmonic components. 

2) Application of circumferential harmonic-dependent inter-segment com- 
patibility constraints. 

3) Solution of frequency-dependent circumferential harmonic components of 
displacements. 

4) Recovery of frequency-dependent response (displacements, stresses, loads, 
etc.) in various segments of the total structure. 
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An overall flowchart outlining the solution algorithm is shown in Figure 4. 
Provision to include the differential stiffness due to the steady loads is also 
shown. 

1 . Transformation of Applied Loads 

The transformation to frequency-dependent circumferential harmonic components 
depends on the form in which the excitation is specified by the user. The follow- 
ing options are considered in the present development to specify the form of exci- 
tation due to the directly applied loads and base acceleration loads: 

Directly applied loads specified as: 

- periodic functions of time on various segments, 

- periodic functions of time for various circumferential harmonic indices 

- functions of frequency on various segments 

- functions of frequency for various circumferential harmonic indices. 

Base acceleration specified as: 

- function of frequency for circumferential harmonic indices 0 (axial) 
and 1 (lateral). 

Details of each of the above five loading conditions are as follows: 

Directly applied loads (segment-dependent and periodic in time) 

If represents a general periodic load on sector n specified as a function 
of time at M equally spaced instances of time per period (Figure 5 ), the lead 
at m^^ time instant can be written as 

m -0 -£c -is 

pf' = p" + ); p" cos(nTT£b) + P" sin(mT£b) 

£=1 ^ 

m = 1 , 2 , . . . , M , 

where b = 2v/l], L = (M-1 )/2 for odd M, Z, = (H-2)/2 for even M. The last term 

L u 

in equation (1) exists only when H is even. The coefficients pP 
("t" = 0; £c, £s, t=l, 2, ..., M/2) in equation (1) are independent of time, 

and are defined by the relations 
*0 n 

P" = i y P" , («. = 0) Part of (2) 

m=l 


-M/2 

+ (-1)"'"V" , (1) 
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-Ic 

,n _ 2 

•Is 
,n 2 


I cos(m-1£b), 
tn=l 


m 


'' OF POOR QUALITY I 

2 , .... 


4 I £in(m-Hb), and 
^ m=l 

-M/2 „ m 

P" = Pi I (-1)'""^ P" (M even only) (a=M/2). 


> ( 2 Contd. ) 


tn=l 


•1 on 


Each of the coefficient vectors P on the left hand sides of equations 
(2) can further be expanded in a circumferential (truncated) Fourier series 


k. 


p" = p + y 

k=l 0 


«l j. II » " 

P*'^ cos (rrTka) + sin(n-lka) 


where n - 1, 2, .... N , 

= 0; £c, JLS.t =1,2, ..., M/2 

a = 2 tt/N 


(3) 


(4) 


kj_ = (N-D/2 for N odd 
k|^=(N-2)/2 for N even. 

The last term in equation (3) exists only when N is even. The Fourier 

coefficients P ^ ("k" = 0; kc, ks, k = 1 , 2, . . . , k^^; N/2) in equation (3) 
do not vary from sector to sector, and are defined by 


N 




-"r 

^ = TT I P" 

n=l 


-" 5 ." 




(k = 0) 


P^*^ = -I y P'^ cos(n-lka) 
n=l 


(k = 1, 2. 


... k^) > (5) 




<s 2 
■ N L 


-"t" 


I P*^ sin (n-lka), and 
n=l 


-"S'" 


^ I P*^ (‘^ even only) (k - N/2) . 


n=l 
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* * t. f. • » • 9 


J'o*' 
— "k*' 

n ^ 


The terms P ("o" = 0; 2 ,c, »s, £ = 1, 2, .... £j_; M/2 and "k" = 0; kc, ks, 
k = 1, 2, .... kj^; N/2) are the transformed frequency-dependent circumferential 
harmonic components of the directly applied loads pR (m = 1, 2, M and n - 1, 

2, .... N). 


Directly applied loads (Cir c umferential harmonic-dependent and periodic in 


time) . 


Such loads can be represented as 


m -0 I, 
P = p ^ + I 


t-1 L 


-£c 

_n I. II 

n 


-£S 

r> ^ . 


P cos(m-Hb) + p sin (m-l£b) 


-H/2 

+ (-1)"’“''/'''" , (6) 


where m = 1, 2, .... M represent the time instances at which harmonic components 
"k" = 0; kc, ks, k = 1, 2, ..., k, ; N/2 of directly applied loads are specified. 

•The coefficients -^"k" on the right hand side of equation (6) are obtained 
using equations (2) with sector number n replaced by harmonic number "k". 

Directly applied loads (freqiiency-and segment-dependent) 

This type of loads can be represented as 


„n 


= P^ 


k^ 

+ 1 ' 
k=l 


- « f) IS 




cos(n-lka) sin(n-lka) 


(-l)n-l ^/2 


(7) 


where "t" (=1, 2, ..., F) now represents the frequencies at which excitation is 

specified. The transformed frequency-dependent circumferential harmonic components 
-"t" 

M|.M • 

P ("k" = 0; kc, ks, k = 1 , 2, ..., k^^; N/2) are obtained using equations (5) 
with as defined above. 


Directly applied loads (frequency- and circumferential harmonic-dependent) 

These loads are the transformed frequency-dependent circumferential harmonic 
components ("k" = 0; kc, ks, k = 1, 2, kj^; N/2) with (=1, 2, ..., F) 

representing the various frequencies at which the directly applied loads are 
specified. 

Base acceleration (frequen cy- and circumferential harironic-dopendent) 

In Appendix A, it is shown that the components of the translational base 
acceleration contribute to inertial loads on the rotating structure in the follow- 
ing manner: 

l,. i ■— 

CF POOR QUALITY' 
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1. Axial component contributes to P where "k" = 0, and represents the 
specified excitation frequencies. 

HI. II 

2. Lateral components contribute to P where "k" = Ic and Is, and "Z“ re- 
presents the effective excitation frequencies which are shifted from the 
specified frequencies by + the rotational frequency. 

2. Application of Inter-Segment Compatibility Constraints 


As shown in Section 4.5.1 of Reference 7, equations (2) of Section 2 
are used to derive the compatibility conditions relating the circumferential 
harmonic component degrees of freedom on the two sides of a rotationally cyclic 
sector: 


side 2 side 1 


=1 

■ =? 




= u^^ cos(ka) 


sin(ka) 


= -u^'" sin(ka) 


cos(ka) 


= 

Ui . 




(k = o) 


S (k = 1, 2, .... k^) > (8) 


(k = H/2) 


In order to apply these constraint relationships for any given harmonic k, 
an independent set u consisting of the circumferential harmonic component (cosine 
and sine) degrees of freedom from the interior and side 1 of the cyclic sector is 
defined, u is selected from the 'analysis' set degrees of freedom (i.e., the 
degrees of freedom retained after the application of constraints and any other 
reduction procedures), and is defined as 


= Gsk(k) 


> 


(9) 


u^*' and u*^^ each contain all (and only) the 'analysis' set degrees of freedom 
from the interior and both sides of the cyclic sector. Equations (8) are used 
to define some of the elements of the transformation matrices «ck and G^j,. For 
k = 0 and N/2, the matrix is null. 


OF POOH QUALITY 
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3. Solution of Frequenc'>Bf'-?!nd5nt Hr'rnonl Displacements 


'< c ruiL S'i'4S;>r* 

OF POOR QUALITY 


For a given harmonic k, the introduction of IT' in the equations of motion, 
equations (38), Section 2 , results in the transformed equations of motion (Ref. 1 ) 

, ( 10 ) 
where # = h" t h" G^,. , 

^ =ck ♦ 4 Sk • 

L (H) 

' =Jk =ck " <k >=" Ssk • 

t 

As discussed in subsection 1 of Section 3 , P and P are the transformed 
frequency-dependent circumferential harmonic components of the directly applied and 
base acceleration loads. 

At any excitation frequency o)", let 


^ and accordingly, 

“K =K ^ 

where P and u are complex quantities. Equation (10) can be rewritten as 


r J2. , •r-K-i-K gK 

L-tO M + loj 6 + K JU = P 

The excitation frequency tc' is given by 


w' = to for all directly applied and axial base acceleration 

loads, and . (14) 

= to+n for lateral base acceleration loads. 

Equation (13) is solved for u for all excitation frequencies and all harmonics 
as specified by the user. The cosine and sine harmonic components of displacements 
are recovered using equations (9) . 

4 . Recovery of Frequency-Dependent Displacements in Various Segments 

This step is carried out only when the applied loads are specified on the 
various segments of the complete structure. 


- 23 - 



For loads specified es functions of time, equation (3) is used to obtain the 

•••’ ^max* 
The user specifies 


displacements in various uognsents with "t" = 0; tc, is, i. 


The circumferential harmonic k is varied from k^. to k 

mm max 

and 

max min max 

For loads specified as^functions of frequency, equation (7) is used to 

obtain the displacements n" in various segments with representing the 

u 

excitation frequencies. Tlie circumferential harmonic is varied from user speci- 

fied k„. to k„,„. 
mm max 

The solution procedure discussed in this section has been implemented as a 
new capability in MASTR.AH (Ref. 8 ). 
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4. EXAJ1PLES 


Five inter-related examples are presented to illustrate the theoretical 
development of the previous sections. The nevi capability added to NASTRAN to 
conduct forced vibration analysis of rotating cyclic structures (Ref. 8 ), 
as a result of the present development, has been used to conduct these examples. 

A 12-bladed disc is used for demonstration. 

Example 1 is conducted on a finite element model of the complete structure 
(Figure 6 )- Examples 2 through 5 use a finite element model of one rotationally 
cyclic sector (Figure 7 ). Results of example 1 are used to verify some of the 
results obtained in the remaining examples. Table 1 summarizes the principal 
features demonstrated by these examples. 

Steady-state frequency-dependent (sinusoidal) or time-dependent (periodic) 
loads are applied to selected grid point degrees of freedom. The specified loads 
can represent either the physical loads on various segments or their circumfer- 
ential harmonic components. For illustration purposes only, the frequency band 
of excitation, 1700-1920 Hz, due to directly applied loads and base acceleration 
is selected to include the second bending mode of the disc for a circumferential 
harmonic index k = 2. The 'blade-to-blade' distribution of the directly applied 
loads also corresponds to k = 2. Table 2 lists the first few natural frequencies 
of the bladed disc for k = 0, 1 and 2. Modes for k = 2 are shown in Figure 8 . 


General Input 


1 . Parameters : 


Diameter at blade tip = 
Diameter at blade root = 
Shaft diameter = 
Disc thickness = 
Blade thickness = 
Young's modulus = 
Poisson’s ratio = 
Material density = 
Uniform structural damping 


19.4 in. 

14.2 in. 

4.0 in. 

0.25 in. 

0.125 in. 

30.0 X 10® Ibf/in^ 

0.3 

7.4 X 10'^ Ibf-sec^/in^ 

(g) = 0.02 


2. Constraints: 

All constraints are applied in body-fixed global coordinate system(s). 
All grid points on the shaft diameter are completely fixed. Rotational degrees 
of freedom 9^ at remaining grid points are constrained to zero. 
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Figure 6: NASTRAN tiodel of the 12-Bladed Disc 
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TABLE 1: PRINCIPAL FEATURES DEMONSTRATED BY EXAf^PLE PROBLEMS 


Example 

No. 

Finite 

Element 

Model 

of 

Applied loads specified as functions of 

Base 

Acceleration 

Rotational 

Speed 

Frequency (sinusoidal) 

Time (periodic) 

Physical 

Components 

Circum. Harmonic 
Components 

Physical 

Components 

Circum. Harmonic 
Components 

1 

Complete 

Structure 

Yes 




No 

No 

2 

Cycl ic 
Sector 

Yes 




No 

No 

3 

Cyclic 

Sector 


Yes 



Yes 

Yes 

4 

Cycl ic 
Sector 



Yes 


No 

Yes 

5 

Cycl ic 
Sector 




Yes 

No 

Yes 


O O 
•n :?3 

T3 O 

8 

SJ {••' 

<o n 

cz : -j 
> {. • 
n i.; 

•H Tt 

•< C;> 



_ i 




a I » .C..« 

TABLE 2: DLAOED-DISC fiATURAL FREQUENCIES POOS QUALITY 


Frequency (Mode No.), 

Hz. 

Mode Description 

. 

o 

li 

k = 1 

■a 

214 (1) 

208 (1) 

242 (1) 

T 

591 (2) 

594 (2) 

622 (2) 

1 

1 aaBBCtfr- ■ 

1 

1577 (3) 

1633 (3) 

1814 (3) 

1 

2468 (5)* ** 

2460 (4) 

2433 (4) 

! ^ 


* k is the circumferential harmonic index 

** Mode No. 4 for k = 0 at 1994 Hz represents an in-plane shear mode not excited 
by the applied forces. 
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EXAtiPLE 1 • 

Description 

This example uses the direct frequency response capability in NASTRAN, 

RF8, and forms the basis to verify some of the results of examples 2 through 5. 

Input 

1. Parameters: 

Same as general input parameters. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

P(f;n) = A(f) cos ) . 

where n is the segment number, 

(D represents k = 2, 

© represents the total number of segments in the bladed 
disc. 

P is specified using RLOADi bulk data cards. 

Results 

Sample plots of grid point displacement and element stress response are 
shown in Figures 9 through 11. The expected behavior about a k = 2 natural 
frequency of the bladed disc can be seen in all these figures. 

All the response plots (Figures 9 through 25 except 16) have been obtained 
using the plotting capability of NASTRAN. On any given plot, the various curves 
are identified, in order, by symbols X, *, +, . and o. The sequence of 

curves is indicated in the first line of the plot description at the bottom 
left of each figure. 
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Figure 9 
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Figure 10 
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Description 

This example uses the forced vibration capability with cyclic syrrenatry. 

The user input/output data for loads, displacements, stresses, etc., pertain 
to the physical representation of the various segments of the bleded disc. 

The frequency-dependent applied loads correspond to k = 2, and hence the solution 
loops on the circumferential harmonic index k are restricted to k = 2 only via 
parameters KMIN and W-iAX. 

Input 

1. Parameters; 

In addition to general input parameters, 

CYCIO = +1 physical cyclic input/output data 
KMIN = 2 minimum circumferential harmonic index 
KMAX = 2 maximum circumferential harmonic index 
NSEGS •= 12 number of rotational ly cyclic segments 
RPS - 0.0 rotational speed 

6KAD = FREQRESP\ Specify the form in which the damping parameters 
LGKAD = +1 / are used. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

p"(f) = A(f) cos (7PT-(D.||) . 

where n is the segment number, 

(?) represents k = 2, 

(Q represents the total number of segments in the bladed 
disc. 

P is specified using RLOADi bulk data cards. 

Results 

Displacement and stress output results for selected grid points and elements 
are presented in Figures IZ through ig. Agreement between results of Figures 12-13 
and Figure 9, Figure 1^ and Figure and Figure 15 and Figure 11 is excellent. 
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EXW^PLE 3 

Description 

This example uses the forced vibration capability with cyclic syrrenetry. 

The user input/output data pertain to harmonic representation. Frequency- 
dependent excitation is provided by both directly applied and base acceleration 
loads. 

Input 

1. Parameters: 

In addition to general input parameters, 

CYCIO = -1 harmonic cyclic input/output data 
Kf-ilN = 0 minimum circumferential harmonic index 

KMAX = 2 maximum circumferential harmonic index 

NSEGS = 12 number of rotational ly cyclic sectors 
RPS = 600.0 revolutions per second 

BXTID, BYTID, BZTID \ Refer to TABLEDi bulk data cards to specify 
BXPTID, BYPTID, BZPTIO/ magnitude and phase of base acceleration 

components. 

GKAD = FREQRESP\ Specify the form in which damping parameters are 
LGKAD = +1 / used. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

a) = A(f) specified on RLOADi bulk data cards. 

b) Base acceleration as shown in Figure 16 . 

Results 

Results are shown in Figures 17 through 25. 

Figures 17 and 18 present k = 0 results (subcase 1). The excitation consists 
of axial base acceleration and directly applied loads. The selected frequency 
band of excitation, 1700-1920 Hz, lies between the second out-of-plane disc bending 
mode frequency (1577 Hz, k = 0, Table 2) and the first in-plane shear mode 
frequency (1994 Hz, k " 0, Table 2). Since the excitation is parallel to the 
axis of rotation, only the former mode responds. 
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Figures 19 through 23 present Ic - 1 results (subcases 2 (k = Ic) and 3 
(k = Is)). The excitation is due to lateral base acceleration only. Although 
the frequency band of input base acceleration is 1700-1920 Hz, the rotation of 
the bladed disc at 600 Hz (parameter RPS) splits the input bandwidth into two 
effective bandwidths: 

(1700 - 600) = 1100 to (1920 - 600) = 1^ Hz, and 

(1700 + 600) = 2m to (1920 + 600) = 2520 Hz. 

The only k = 1 mode in these effective band-widths is the first torsional 
mode of the blade vfith the disc practically stationary (2460 Hz, k = 1, Table 2). 

This is shown by the out-of-plane displacement magnitudes of grid points IS 
(blade) and 8 (disc) respectively (Figures 19 (k = Ic) and 22 (k = Is)). The 
corresponding phase responses of these grid points are sho-wn in Figure 21 . 

Figures 24 and 25 present k = 2 results (subcase 4 (k = 2c)). The excitation 

consists of directly applied k = 2c loads. The out-of-plane displacement magni- 
tude of grid point 18 (Figure 24) compares v<ell with that obtained in example 2 
(Figure 13 . Table 3 lists the out-of-plane displacement response of grid point 
18 as obtained in examples 2 and 3. The marginal difference in response in example 3 
is due to the Coriolis and centripetal acceleration effects at a rotational speed 
of 600 revolutions per second. 

No k = 2s loads are applied in this example (subcase 5). 
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TABLE 3; EFFECT OF CORIOLIS AtJD CENTRIPETAL ACCELERATIONS ON THE 
DISPLACEMENT RESPONSE OF GRID POINT 18 AT 600 RPS. 


Frequency 

Hz 

Example 2 

Example 3 I 

Segment 1 (subcase 1 ) 
Mao. (in)/Phase (deq) 


1700 

7.2655 E-5/349.4 

7.5132 E-5/354.3 

. 

1750 

1.3071 E-4/343.1 

1.3844 E-4/347.3 

1770 

2.1580 E-4/332.7 

2.3252 E-4/335.8 

1796 

3.4139 E-4/314.6 

3.7252 E-4/315.2 

1814 

4.8374 E-4/269.9 

4.9177 E-4/266.8 

1832 

3.4146 E-4/224.9 

3.2655 E-4/225.5 

1850 

2.1451 E-4/206.6 

2.0742 E-4/209.3 

1880 

.1.2433 E-4/195.6 

1.2214 E-4/199.2 

1920 

7.6125 E-5/190.4 

7.5397 E-5/194.3 




\ 
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EXAMPLE 4 


Description 

This example uses the forced vibration capability with cyclic syniiietry. 

The user input/output pertains to physical representation, periodic loads are 
specified as functions of time on the segments of the bladed disc corresponding 
to k = 2. For clarity of illustration only, sinusoidal loads of varying 
amplitudes at a frequency of 1814 Hz are specified. The Fourier decomposition 
of these sine functions obviously contains contributions from first harmonic 
alone (£ = 1)-- the parameter LMAX accordingly has been set at 1 (£. = 0, Ic, Is). 

Input 

1. Parameters: 

In addition to general input parameters, 

CYCIO = +1 physical cyclic input/output data 

KMIN = 2 minimum circumferential harmonic index 

KMAX = 2 maximum circumferential harmonic index 

LMAX = 1 maximum harmonic in the Fourier decomposition of periodic, 

time-dependent loads, 

NSEGS = 12 number of rotational ly cyclic sectors 
RPS = 600.0 revolutions per second 

GKAD = FREQRESPj Specify the form in which the damping parameters are 
LGKAD = +1 / used. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

P"(t) = A(t) cos (]rr-(2>||) , 

where n is the segment number, 

(^represents k = 2, 

(Jj) represents the total number of segments in the bladed 
disc, 

A(t) = A*sin (2TT*1814*t). 

P is specified on TLOADi bulk data cords. 
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T 

Results 

Results are presented In Table 4 and are in good agreement with those from 
example 3. 

r . 






TABLE 4; COMPARISON OF RESPONSE AT 1814 Hz 


Grid Pt.Disp. or 
Elem. Stresses 

1 E>;amDle 3 

Examole 4 

Example 5 


Segment 1 (subcase 1) 
Maq. ( in)/Ph3se(dsq) 

k - 2c'Xsubcase 4) 
Maq.(in)/Phase(deq) 

8 (T3RM), iJ^ 

5.4297 E-4/82.6 

5.4299 E-4/82.6 

5.4299 E-4/82.6 

18 (T3R?D, 11 

4.9177 E-4/266.8 

4.9180 E-4/266.8 

4.9180 E-4/266.8 

" <3>- V,!* 

1.4841 E 3/84.7 

1.4842 E 3/84.7 

1.4342 E3/84.7 

n (5). 

2.0891 E 2/83.4 

2.0892 E 2/83.4 

2.0892 E2/83.4 

” v.i 

1.0774 E 2/64.7 

1.0775 E 2/64.7 

1.0775 E2/64.7 

" <"»■ “xx.a* 

1.4677 E 3/263.3 

1.4678 E3/263.3 

1.4678 E3/263.3 

n oa). 

2.2489 E 2/260.3 

2.2491 E 2/260.4 

2.2491 E2/260.4 

11 (14). 

1.8510 E 2/253.0 

1.8511 E 2/253.0 

1.8512 E2/253.0 


* Fibre distances 1 and 2. 
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Description 

This example uses the forced vibration capability vn’th cyclic syrtraetry. 

The user input/output pertains to hanaonlc representation. Periodic loads are 
specified as functions of time for the circumferential harironic index k = 2. 

For clarity of illustration only, sinusoidal loads are selected. 

Input 

1. Parameters: 

In addition to general input parameters, 

CYCIO -1 harmonic cyclic input/output data 
KHIH = 2 minimum circumferential harmonic index 
KJ1AX = 2 maximum circumferential harmonic index 
LMAX = 1 maximum harmcnic in the Fourier decomposition of periodic, 
time-dependent loads. 

NSEGS = 12 nurber of rotaticnally cyclic sectors 
RPS = 600.0 revolutions per second 

6KAD = FREQRESP\ Specify the form in which the damping parameters 
LGKAD = +1 / are used. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

P^'^(t) = A-sin (2Ti-1014-t) , 

specified on TLOADi bulk data cards. 

Results 

Results are presented in Table 4 and agree well with those from example 3. 
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5. COiJCLUSlOHS 

1. A nevj capability has boon dovelcped and added to the general purpose 
finite element program NASTfWfl Level 17.7 to conduct forced vibration analysis 
of tuned cyclic structures rotating about their axis of symmetry. 

2. The effects of Coriolis and centripetal accelerations together with 
those due to the translational acceleration of the axis of rotation have been 
included. 

3. A variety cf user options is provided to specify the loads on the 
rotating structure. 

4. Five interrelated examples are presented to illustrate the ,'arious 
features of this development. 
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S. RECCM’IEKDATIONS 


1. This IS a new capability and, therefore, the examples presented herein 
have been primarily designed to illustrate the various basic features of the 
development. Application to a variety of real problems would substantially 
contribute twards determining its merits and limitations with regards to its 
applicability, usefulness, and savings in modelling and computational time. 

2. The capability should be extended to conduct forced response analysis 
using noriral modes with cyclic symmetry as the basis. 

3. Inclusion of induced and applied oscillatory aerodynamic loads within 
the capability would be a desirable step in solving the forced vibration problems 
of turbomachines. 
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APPENDIX 

INERTIAL LOADS DUE TO BASE ACCELERATION 

The acceleration of the axis of rotation generates inertial loads at all 
grid points of the complete structure. In this appendix, the generation of these 
inertial loads and their transformation to frequency-dependent circumferential 
harmonic components are discussed. 

As given by equation (27) of Section 2 , the inertial forces on the 
three translational degrees of freedom at an arbitrary point P of the modelled 
cyclic sector, expressed in the global (displacement) coordinate system, are 

{P®} = -EH23CRo^ " ’ (1) 

where 
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with c H cos fit and s = sin fit. 

Since all the cyclic sectors are identical in all respects except for the 
specified loads, no generality is lost in assuming, for simplicity, that the 
modelled sector is the n = 1 sector. Equation (1) can, then, be rewritten as 

{P^} (3) 


where 

c^ = cos (n'-l • 1 • 27r/N), and 

s^ = sin (rf-1 • 1 • 2r/N) 



{pG) = [t^Gj 


0 0 


-s. 


'n J 


Substituting equation (3) in equations (5) of Section 


and noting that 





I • cos (rirr • k • 2v/u) = N/2 , k=l 


= 0 , kj!1 . 


• 005(71^ • k • 2Tj/rj) = 0 


(5) 

(coritd) 


• sin(rPT • k • Zir/tl) = 0 , 


' • sin(rTT • k • 2 .t/H) = N/2 , k=l 

n " 

= 0 . kf[ , 

the circumferential harmonic components of the base acceleration loads become 


{P®} = j ° [ ("k" = 0) 

Lo J . 


{p.C} ^ P 


("k" = Ic) 


{plS} ^ [T®^J ) P 


("k" = Is) , and 


{pkc.ksj ^ ^ 


In the present development, the components of t se acceleration Xg, Yg and 
Zg are considered to be sinusoidal of frequency o', and are specified as 


^0 ^0, mag ’ 


Y = Y c 

’O u, mag 


cos(i‘’t , ar 


h mg ‘ * V ■ 


-F,n- 



From equation (2), therefore, u-a can write 


Py® = ■"’i^^O.mag *J'Y^ ^O.mag ‘ ^ 

= ->"i:-Yo,mag ^ '^yV'’ ^O.rnag J 

The cosine and sine products in equations (8) can be expressed in terms of 
individual cosine and sine terms v/ith frequencies (oo + n) and (w - n). 

The following conclusions about base acceleration loads can, therefore, be 
drawn by substituting equations (8) into equations (6): 

1. The axial component of base acceleration, Xq(w), contributes to ^ at exci- 
tation frequencies to. 

2. The lateral components of base acceleration, YQ(to) and Zg(to), contribute to 
and at excitation frequencies (to ± il) for each to specified. 
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Damping matrix 

Coriolis acceleration coefficient matrix 

Rayleigh's dissipation function 

"Synrnetric Components" transformation matrix 

Unit vectors along Inertial XYZ axes 

Unit vectors along Basic Xg Ig Zg a>^es r (Figure 3) 

Unit vectors along Global xyz axes J 

Stiffness matrix 
Circumferential harmonic index 
Tima harmonic index 

Hass matrix, number of time intervals per period (Figure 5) 
Centripetal acceleration coefficient matrix 
Base acceleration coefficient matrix 
Hass 

Number of cyclic sectors in the complete structure 
Load vector 

Aerodynamic coefficient matrix 
Base acceleration vector 
Position vectors (Figure 3) 

Kinetic energy, coordinate systen tranofortnation matrix 
Time 

Strain energy 

Physical displacement degrees of freedom 

Virtual work 
Rotational frequency 
Forcing frequency 


I 
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SVi-^BOLS {Continued) 


Superscripts 


B 

G 

K" 


m 


n 


-0 


'k 


-zc 


-zs 

-kc *• 
-ks 
-H/2 
-N/2 


Basic 

Global 

Independent solution set in "syrmietric components" 
nth time instant 
nth cyclic sector 


Fourier coefficients ("symmetric components") 
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